Sliding Phases in XY-Models, Crystals, and Cationic Lipid-DNA Complexes 
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We predict the existence of a totally new class of phases in weakly-coupled, three-dimensional stacks 
of two-dimensional (2D) XF-models. These "sliding phases" behave essentially like decoupled, 
independent 2D XY-models with precisely zero free energy cost associated with rotating spins in 
one layer relative to those in neighboring layers. As a result, the two-point spin correlation function 
decays algebraically with in-plane separation. Our results, which contradict past studies because we 
include higher-gradient couplings between layers, also apply to crystals and may explain recently 
observed behavior in cationic lipid-DNA complexes. 
Pacs: 61.30.Cz,61.30.Jf,64. 70.Md 
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Spatial dimensionality greatly affects the nature of or- 
der in condensed matter systems. Three-dimensional 
(3D) ATY-systems such as superfluids and ferromagnets 
have true long-range order with divergent correlation 
lengths at a second-order transition separating the high- 
temperature disordered phase from the low-temperature 
ordered phase. Two-dimensional (2D) ATY-systems, in 
contrast, exhibit power-law decay of correlations in the 
low-temperature phase. At high temperatures beyond 
the Kosterlitz-Thouless (KT) transition temperature, 
thermally excited vortices destroy the quasi-long-range 
order and cause correlations to decay exponentially . 

Many experimentally realizable systems such as lay- 
ered superconductors [0], free-standing liquid-crystal 
films [0] , and lyotropic smectics with internal membrane 
order can be viewed as stacks of two-dimensional layers 
with interlayer couplings Q that can be varied substan- 
tially, for example by changing the layer spacing. What 
is the phase behavior of such a system? If there is no 
coupling between layers, each will exhibit 2D behavior; if 
the coupling is strong, the system will exhibit 3D behav- 
ior. Shortly after the discovery of the KT transition, it 
was suggested that a weakly-coupled stack of XY-models 
might behave over some temperature range as a stack of 
decoupled layers, i.e., that such a system could proceed 
from three-dimensional behavior at low temperatures to 
a 2D power-law phase at intermediate temperatures to a 
disordered phase at high temperatures . 

Subsequent studies, however, demonstrated that if the 
only interlayer couplings are Josephson (i.e., proportional 
to cos(6>„ — 6 n +i), where n is the A"Y-angle variable 
in layer n), the intermediate 2D power-law phase is 
squeezed out, and the system goes directly from the 3D 
long-range ordered phase to the disordered one with in- 
creasing temperature. This happens because the "decou- 
pling temperature" Td above which the Josephson cou- 
pling becomes irrelevant is greater than the Kosterlitz- 
Thouless temperature Tkt below which the 2D ordered 



phase is stable against vortex unbinding. Thus, the tem- 
perature window Td < T < Tkt over which the 2D slid- 
ing phase can exist disappears. 

In this paper, we revisit this old and seemingly dead 
issue and show that a thermodynamically stable phase 
exhibiting 2D power-law correlations is in fact possible. 
The new ingredient in our analysis, which was not present 
in previous treatments, is competing higher-order gradi- 
ent couplings between layers |7j. These gradient cou- 
plings, in the absence of Josephson couplings between 
layers, produce two-point correlation functions that are 
identical in form to those of a stack of decoupled 2D 
layers. We will refer to this phase as a sliding and not 
a decoupled phase because the ATY-angle variables in 
different layers can slide relative to each other without 
changing the energy of the system and because nonzero 
couplings between layers (though not of the Josephson 
type) are, in fact, present in our model, and furthermore, 
necessary for the existence of this phase. Remarkably, it 
is possible through judicious tuning of interlayer gradi- 
ent couplings to satisfy Td < Tkt and produce a stable 
sliding phase for T d <T < T K t- 

This investigation into whether or not a sliding phase 
of XY-models exists was inspired by recent work on the 
possible sliding columnar phase in cationic lipid-DNA 
complexes ||. In these complexes, DNA molecules are 
intercalated between lipid bilayers and, within each layer, 
the molecules are situated on a one-dimensional lattice. 
Experiments may be consistent with the existence of a 
sliding columnar phase in which lattices in neighboring 
layers are able to slide over each other without energy 
cost || in complete analogy to the sliding phase just de- 
scribed for XY-models. Indeed, we have shown theoreti- 
cally, by methods analogous to those presented here, that 
it is possible to have a sliding columnar phase in these 
complexes and, in addition, possible to have a sliding 
phase in a layered crystal. Details on these two systems 
will be presented in a future publication |lO|; for the re- 
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mainder of this paper, we will focus on XF-models. 

The traditional theory for a stack of XY- models begins 
with sum of independent Xy-Hamiltonians 



K x 



(1) 



where r = (x,y,0) is a point in the x-y plane and Vj_ 
is the gradient operator acting on these two coordinates. 
Josephson-like couplings between layers are then added; 
these are given by 



Hj[s n ] = -Vj[s n ] / d 



V 



s P 9 n +p(r) 



(2) 



where s n is an integer-valued function of layer number 
n satisfying ^ n s n = if there arc no external fields 
inducing long-range order. If all Vj[s„] are zero, then 
in the low-temperature phase (6^(r))o = rj\og(L/b) and 
cos[0„(r) - n (O)] ~ r^, where 



V 
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(3) 



L is the sample width, b is a short-distance cutoff in 
the x-y plane, and (.)q refers to an average with re- 
spect to Ho . The averages of the Josephson Hamiltonians 
with respect to 7i scale as (Hj[s n ]) ~ £ 2 -'>[ s ™l where 
V[ s n] = 2^2p s i>- Clearly, the most relevant Josephson 
coupling is the one with the smallest value of 7?[s n ], which 
results when s n is non-zero on the smallest number of 
planes. Since J2 n Sn = 0, the smallest value of T)[s n ] is 
obtained for couplings between two layers separated by 
p layers with s n = = $n.o ~ $n,p- For these two-layer 



couplings, rf 



77 [s£ 



77 for all values of p. Thus, 



the decoupling temperature above which all Josephson 
couplings are irrelevant is Td = AirK. The 2D KT tran- 
sition for decoupled layers is Tkt = ifK/2\ this implies 
Tkt/Tci = 1/8 < 1, and there is no decoupled phase with 
power-law correlations. 

Josephson couplings are not, however, the only ones 
permitted by symmetry. Gradients of 9 n in different lay- 
ers may also be coupled. The Hamiltonian for the ideal 
sliding phase is Hs = Ho + H g , where 



H a 



\Y,f*r?? ~ 0n(r))} 2 . (4) 



This Hamiltonian is invariant with respect to 9 n (r) — > 
9 n (r) + ip n for any constant ip n , i.e., the energy is un- 
changed when angles in different layers slide relative to 
one another by arbitrary amounts. The sliding Hamilto- 
nian can be written as 



(•5) 



where K nn > = Kf n - n , with /„ = (1 + Z) m 7ro)<W) - 
) and 7 m = U m /K. Also, the 

Fourier transform 



f(k) = 1 + ^2 7 m (l - cos ki 



(6) 



of the reduced coupling /„ will be used extensively below. 

Correlations in the sliding phase can easily be calcu- 
lated from Eq. (J5J) - We find 

<0„(r)Mr')> s = Vf-lJHL/b) E(\r r'|)], (7) 

where (-)s is an average with respect to Hs and E(r) = 
f dq[l — J (qr)]/q tends to zero as r — > and to ln(r/6') 
with b' /b m 0.2 as r — > 00. The inverse coupling ff 1 is 
defined by 



/p ~W f(k) ■ 



(8) 



Thus we find (0*(r)) s = 77^(0) ln(L/6) and 

gs(r,p) = ([9 n+ p(r)-0 n (O)} 2 ) s 

= 2[?j s (p) HL/b) + Vs (p) Hr/b')} (9) 

for large r. The coefficients of the logarithms are 

Vs(p) = Vifo 1 - fp 1 ) and Vs (p) = riff 1 . (10) 

Note that 77s (p) — r}5p,o and rjsip) = 7/(1 — 5 Pi o) when 
H g = 0. Using Eq. we find that the correlation func- 
tion Gs(r,p) = (cos[0„ +p (r) — 0„(O)]) S satisfies 



Gs(r,p) 



{ (L/by 

\ (r/V) 



-Vs(0) 



P^Q 
p = 0. 



(11) 



Thus, the two-point spin correlation function for spins 
in different layers vanishes in the L — > 00 limit, whereas 
that for spins in the same layer has exactly the same form 
as it would for a stack of decoupled layers. Now, how- 
ever, the exponent 775(0) depends on the detailed form of 
the interlayer gradient couplings via f(k). The two-point 
spin correlation function is zero for spins in different lay- 
ers; nonetheless, nonvanishing couplings between layers 
cause other correlation functions that are zero for the 
totally decoupled layers to become nonzero. 

Having established that the sliding phase (if it has 
not melted) behaves like a stack of decoupled layers, 
we now ask what happens when the Josephson inter- 
layer couplings of Eq. (^) are turned on. From Eq. 
(@), we find that {Hj[s n ]) s ~ L 2 "^[^l where ?} s [s n } = 
n ' s n s n' ffl n i ■ As for the decoupled case, the mini- 
mum value of rjs [s n ] is obtained when s n = s^. Thus the 
decoupling temperature for couplings with s n = is 



Tdip) 



fo fp 



(12) 



which depends on p. The temperature above which all 
Josephson couplings are irrelevant is = max p T^(p). 
We will show later that this maximum over all p is finite. 
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To prove the stability of the sliding phase, we must 
show that, for some range of the couplings U m , the de- 
coupling temperature T d calculated above is less than 
Tkti the temperature at which vortices unbind. To cal- 
culate Tkti we must calculate the vortex energy. This 
calculation in our model is similar to that for decoupled 
layers. Vortex excitations in individual layers remain well 
defined when the layers are coupled, although, when cou- 
plings are sufficiently strong the system becomes truly 
three dimensional, and vortices should be viewed as seg- 
ments of closed vortex loops. Defining v„ (r) — Vj_0 n (r), 
we have 



* v„ • dl = 2n V" k n j, 
Jr , 



(13) 



where k n> i is the integer strength of the Ith vortex in 
the nth layer and T is a contour in layer n enclosing the 
vortices. Applying Stokes theorem to Eq. ( |l3| ) then gives 



where 



V_l x v„ = m z n (r)z, 



<(r) = 2TrJ2k nA 6 2 (r-r ntl ) 



(14) 



(15) 



is the vortex density in layer n and r nA gives the position 
of each vortex in the x^y plane. We then take the 2D curl 
of both sides of Eq. (|14| ) and Fourier transform to find 



^n(q_L) 



" >r<l -j >"',.'■ <\ 



(16) 



Using this result in Eq. (|^), we obtain the vortex energy 
d 2 q± _ m^(q_L)m^(-q_L) 



E 



9 /"- ? 



(2n)< 



= nK /»-»' ( J2 fc «>< ) ( E W ) 
n,n' \ I / \ V ) 

-nK ^2 fn-n'k n ,ik n ',i'E(\r nt i -r n > t i>\). (17) 

Since /„_„' is a positive definite matrix, this equation 
implies that in the thermodynamic limit, there must be 
charge neutrality in each layer, i.e., k n j = 0. The 
interaction between like-sign vortices in different layers n 
and n' is attractive if f n -n' < and repulsive if f n -n' > 
0. Since we assume that individual layers are stable in the 
absence of couplings between layers, /o > and like-sign 
vortices within a single layer repel. 

The vortex energy Ey and Boltzmann statistics im- 
ply that the number of times a given configuration of 

vortices occurs in the system scales with system size as 
£2-wtK] where 



(18) 



a„ = J2i and the factor of L 2 counts the number of 
places in the 2D plane the configuration can be placed. 
Clearly, if r)[a n ] < 2, the particular vortex configuration 
{<j n } will proliferate. The "Kosterlitz-Thouless unbind- 
ing" for {cr„} therefore occurs at a temperature 



ttK 

n.n 1 



fn—n' &n &w 



(19) 



If there is only one vortex in layer 0, then a n = <r„ = 5 Jly o 
and TktIo'u] = TrKfo/2. If there is a +1 vortex in layer 
zero and ail vortex in layer p, a n = cr£ = S n fi ± S n . p 
and TktW?^] = n K (fo ± fp)- Note that when f p is 
nonzero TktI^^ 1 ] is not twice Tkt[<^^\- In fact it is 
possible for Tkt^h] to be less than TktWu] f° r one or 
more configurations {cr n }. The interactions between lay- 
ers lead to composite multi-layer vortices that cost less 
energy to create than a single vortex in an individual 
layer. Unbinding of bound pairs of any set of individual- 
layer or composite vortices will destroy the rigidity within 
those layers. Thus, the transition temperature to the dis- 
ordered state is Tkt = m hi{ CTre } Tkt\gt\, and the sliding 
phase exists provided 



T KT min^ T KT [a n ] 



m&x p T d (p) 



> 1. 



(20) 



We will now discuss how the interlayer gradient poten- 
tials U m can be chosen so that (3 > 1. The basic strategy 
is to choose the U m so that f{k) has a minimum near 
zero at some value of k. We consider a model with both 
first- and second-neighbor interactions. We ensure that 
there is a minimum in f(k) at kg by requiring 

/(fc ) = 1 + 7i (1 - cos ko) + 72 (1 - cos 2fc ) = A (21) 

and /'(fco) = 0. These two conditions determine 71 and 
72 in terms of ko and A. In the range of ko and A we 
consider, 71 > > 72 and 71 > I72I. The minimum 
can be tuned to zero by taking A to zero, in which case 
f(k ) = but f(k) > for all k ^ k . For small A, 
Jq 1 — f^ 1 is dominated by values of k near ko, and we 
have 



f -i _ f -i 
■10 Jp 



1 - cos(pfc )e' p \ / ^ 



P_ 

C 



VCA 
(pk -2nl) 2 



(22) 



where the final form is valid for pko ~ 2ttI, I is an inte- 
ger, and C = f"(k )/2. From Eq. fl22]), we see that there 
exists a curve T d (p, k ) ~ [p/C + {pk - 2tt/) 2 /2^/AC}- 1 
for each value of p that specifies the decoupling temper- 
ature as a function of ko- For fixed fco, max p T d (j>, ko) 
occurs at p = [27r//fco] if < {2itl/ko\ < 1/2 and at 
p = [27rZ/fco] + 1 if 1/2 < {2nl/k } < 1, where [x] is the 
greatest integer less than or equal to x and {x} = x — [x] 
is the fractional part of x. As a function of ko near 2irl/p, 
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Td(p, ko) reaches a maximum at ko — 2nl/p and decreases 
sharply away from this point. Also, in the range of ko 
and A we have considered, we can prove that composite 
like-sign vortices in nearest-neighbor planes p and p + 1 
with er„ = 5 n>p + 8 ntP+ i are the first to unbind and thus 
T KT = nK{f + A) = nK{l + 7l /2 + <y 2 ). Since T KT 
is a smooth function of ko, we find that (3 = Tpcr/Td 
has sharply-peaked minima at ko = 2nl/p. Direct eval- 
uation of (3 for A = yields f3 > 1 in the range 
0.24 < k a /n < 0.40 as shown in Fig. @. 
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FIG. 1. /3 = Tkt/Tci is plotted versus ka/n. Local minima 
near ko/it = 21 /p are labeled by (l,p). Other possible integer 
pairs either do not fall in the range 0.24 < ko/n < 0.40 or 
yield larger values of /3 than those shown above. 

Transitions out of the sliding phase are of the 
Kosterlitz-Thouless or roughening type. The transition 
to the high-temperature disordered phase at Tkt is con- 
trolled by K and the fugacity y ++ for composite like-sign 
vortices in neighboring layers. The transition to the low- 
temperature 3D ordered phase is controlled by the first 
V p = Vj[sP] to become relevant and by U p . 

As we have seen, the Josephson couplings V p are irrele- 
vant with respect to the sliding phase for Td <T < Tkt- 
If all V p are set to zero, the two-point correlation function 
Gs(r,p) vanishes for p ^ 0. Even though the V p are ir- 
relevant, they are not zero. They give rise to nonzero 
perturbative contributions to Gs(r,p) even when p is 
nonzero. Consider for simplicity the nearest-neighbor 
Josephson model {V p — VS P: i). Then 

G s (T,p)=(^J J d\ 1 ...d\ p e~^-^/\ (23) 

where $(n,...,r p ,r) = ([A6> (0,r 1 ) + A6> 1 (r 1 ,r 2 ) + . . .+ 
A6» p (r p ,r)] 2 ) s and A6» n (r 1 ,r 2 ) = n (ri) - n (r 2 ). Using 
the fact that the transition at is a KT transition, it 
can be shown that Eq. ( |23| ) yields an exponential decay 
of correlations with Gs(0,p) ~ e~~ p ^° . The interlayer 
correlation length £ z diverges as T — > according to 



£ 2 cx 1/(T — Td), i.e. the interlayer correlation length 
exponent is v z = 1. This divergence signals the devel- 
opment of true long-range orientational order in the 3D 
ordered phase below Td- The derivation of this result will 
be presented in a forthcoming publication flio| . 

The ideas presented here can also be applied to a three- 
dimensional stack of two-dimensional crystals pQl. An 
interaction Hamiltonian analogous to 7i g in Eq. (Q) that 
couples gradients of displacements in different layers can 
be introduced. Power-law exponents and dislocation en- 
ergies again depend on these couplings, and a sliding 
crystal phase between a low-temperature crystalline and 
a higher-temperature hexatic phase [[llj is possible. The 
sliding crystal phase is similar to a model once proposed 
for the smectic B phase in liquid crystals Also, in- 
terlayer gradient couplings for the hexatic angle can be 
introduced to produce a sliding hexatic phase. Thus the 
phase sequence 3D crystal — > sliding crystal — > 3D hex- 
atic — ► sliding hexatic — > disordered layers is in principle 
possible in lamellar systems. 
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